Elastic energy for reflect ion- symmetric 

topologies 

A Majumdart JM Robbins^ & M Zyskin^ * 
School of Mathematics 
University of Bristol, University Walk, Bristol BS8 ITW, UK 

and 

Hewlett-Packard Laboratories, 
Filton Road, Stoke Gifford, Bristol BS12 6QZ, UK 

February 4, 2008 



Abstract 

Nematic liquid crystals in a polyhedral domain, a prototype for 
bistable displays, may be described by a unit-vector field subject to 
tangent boundary conditions. Here we consider the case of a rectangu- 
lar prism. For configurations with reflection-symmetric topologies, we 
derive a new lower bound for the one-constant elastic energy. For cer- 
tain topologies, called conformal and anticonformal, the lower bound 
agrees with a previous result. For the remaining topologies, called 
nonconformal, the new bound is an improvement. For nonconformal 
topologies we derive an upper bound, which differs from the lower 
bound by a factor depending only on the aspect ratios of the prism. 
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1 Introduction 



Present-day liquid crystal displays (eg twisted nematic) are based on monos- 
table cells, wherein, in the absence of external fields, the orientations of the 
liquid crystal molecules assume a single (spatially varying) mean configura- 
tion which is effectively transparent to incident polarised light. To produce 
and maintain optical contrast, voltage pulses, which reorient the molecules, 
must be continually applied. There is considerable interest in developing 
bistable cells, which support two (and possibly more) stable liquid crystal 
configurations with contrasting optical properties. In bistable cells, power 
is needed only to switch between the two states. One mechanism for en- 
gendering bistability is the cell geometry [HIS EI; nematic hquid crystals in 
prototype cells with polyhedral geometrical features (eg, ridges, or posts) are 
found to support multiple configurations. 

As a simple model for such systems, we consider the mean local orien- 
tation of a nematic liquid crystal in a polyhedral domain as described by a 
director field n subject to suitable boundary conditions. The situation we 
consider, strong azimuthal anchoring, is described by tangent boundary con- 
ditions. Tangent boundary conditions require that, on a face of the domain, 
n lies tangent to the face, but is otherwise unconstrained. This implies that 
on the edges of the polyhedron, n is parallel to the edges, and therefore is 
necessarily discontinuous at the vertices. We restrict our attention to direc- 
tor fields which are continuous away from the vertices (ie, as continuous as 
possible). In this case we can unambiguously assign an orientation to the 
director (as the domain is simply connected), and regard n as a unit- vector 
field. 

In jlj, we give a complete topological classification of continuous tangent 
unit-vector fields in a convex polyhedron. An extension to the nonconvex 
and periodic cases, along with a general procedure for analysing a large class 
of such classification problems, is given in In [Hj we obtain a lower bound 
for the one-constant energy in terms of certain topological invariants, the 
trapped areas. The case of a rectangular prism is considered in [7j, where we 
also derive an upper bound for the equilibrium (infimum) energy for a large 
family of topologies called refiection-symmetric conformal and anticonformal. 
For these topologies, the ratio of the upper and lower bounds depends only 
on the aspect ratios of the prism. We also show that topologically nontrivial 
behaviour of configurations close to equilibrium may concentrate near the 
edges, or may be smoothly distributed, depending on the aspect ratios. 

In this paper we consider again the case of a rectangular prism, and 
improve and extend the previous results of 0. Specifically, we derive a new 
lower bound for the energy of refiection-symmetric topologies, expressed in 
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terms of different invariants, namely the wrapping numbers. In general, the 
new lower bound is an improvement on the previous one. We also extend 
the analysis to all reflection-symmetric topologies, not just conformal and 
anticonformal ones. 

While liquid crystal applications are a principal motivation for this work, 
the problems are also of intrinsic mathematical interest. Minimizers of the 
one-constant energy may be regarded as harmonic maps from a Euclidean 
polyhedron to the two-sphere 5*^. The study of harmonic maps between Rie- 
mannian manifolds is an extensive field, and connections to problems in liq- 
uid crystals are well known [H]. For manifolds with boundary, the regularity 
of minimisers for the Dirichlet problem for harmonic maps with sufficiently 
smooth (C^) boundary and Dirichlet data are investigated in 0. However, 
less appears to be known about the case of manifolds with Lipschitz bound- 
ary, eg domains with corners, and for natural, eg tangent boundary con- 
ditions. There are recent strong results on the existence, uniqueness and 
regularity of minimisers for the Dirichlet problem for harmonic maps of fixed 
homotopy type between Riemannian polyhedra for target spaces of negative 
curvature |T^. However, it appears to be much more difficult to obtain cor- 
responding results for target spaces of positive curvature, eg 5^, which we 
encounter in liquid crystals problems. 

The paper is organised as follows. The topological classification of tangent 
unit- vector fields in a rectangular prism is reviewed in Section |21 We intro- 
duce the refiection-symmetric topologies, which are characterised by certain 
invariants - the edge signs e, kink numbers k and trapped area Q - associ- 
ated with one of the prism vertices. In Section IHl we derive a lower bound for 
the one-constant elastic energy. This turns out to depend on the absolute 
values of the wrapping numbers (which may be expressed in terms of e, k and 
Q). For certain topologies, called conformal and anticonformal, for which the 
wrapping numbers all have the same sign, the lower bound can be expressed 
in terms of the trapped area alone, and coincides with the result previously 
derived in [61 Ej . Conformal and anticonformal topologies are characterised 
in Sectional where it is shown that these are precisely the topologies which 
have conformal and anticonformal representatives of the type considered in 
[7]. In Section El we introduce representative configurations for nonconformal 
topologies, and derive from them an upper bound for the elastic energy. This 
differs from the lower bound of Section |21 by a factor depending only on the 
aspect ratios. Appendix ^ contains a derivation of a formula for the kink 
numbers of conformal and anticonformal configurations. 
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2 Reflection-symmetric topologies 



Let us briefly recall the results concerning the classification of continuous 
tangent unit-vector fields n on a rectangular prism P. For convenience, we 
let P denote the prism without its vertices, so that n is continuous on P. 
For definiteness, we take the prism to be given by < < L^, with edge 
lengths Lj ordered so that > Ly > L^. At each vertex of P, denoted 

V = {vx,Vy,Vz), we associate to n a set of topological invariants, namely 
the edge signs, kink numbers and trapped area. The edge signs, denoted ej, 
determine the signs of n on the edges at v relative to the coordinate unit 
vectors, ie 

n(x, Vy, f^) = X, < X < L^, (1) 

and similarly for Cy and eJ. (Of course, this designation is redundant; the 
edge signs eJ and eJ at vertices v and w joined by an edge parallel to the 
j-axis are necessarily the same.) 

The integer- valued kink numbers, denoted kj, count the windings of n 
along a path about v on the face normal to j. By convention, the paths are 
taken to be positively oriented with respect to the outward normal through 
the centre of the face. The minimum possible winding (a net rotation of 
±7r/2) is assigned a kink number of zero. Nonzero windings are designated 
positive or negative according to their orientation with respect to the outward 
normal (either j or —j). The continuity of n (in particular, the absence of 
singularities on the surface of P) implies that the kink numbers on each face 
satisfy a sum rule |3]; for example, on one of the faces F normal to z, it turns 
out that 

E (^^' - li-^r^'^i-lp^'^'eX) = 0. (2) 

Analogous rules hold for the other faces. 

The last invariant, the trapped area, denoted fl^, is the oriented area on 
the unit two-sphere S'^ of the image, n(C"^), of a surface, C^, which separates 

V from the other vertices. That is, letting {6, 0) denote polar coordinates on 

= I sin edOA d(j). (3) 
Expressed as an integral over itself, is given by 

= I T>-&dS. (4) 

Here C"^ is the outward-oriented unit normal on {C^ points towards v) 
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and dS is the area element, while the vector field D(r) is given by 

Dj = ^ejkiidkn X din) ■ n. (5) 

That (jHI) and (jSj) are equivalent follows from the fact that 

D ■ = detc/ncv, (6) 

where dnc^ denotes the Jacobian of the restricted map ncv : ^ S"^ (and, 
as above, is oriented with respect to the outward normal). 

For a rectangular prism, the trapped areas are necessarily odd multiples 
of 7r/2 (the area of a right spherical triangle), and for given values of the 
edge signs and kink numbers, the allowed values of the trapped areas differ 
by multiples of 4tt (whole coverings of the sphere) - see (|T3|l below. The 
continuity of n (the absence of singularities inside P) implies the sum rule 

E = 0- (7) 

V 

One can show ([4 ) that the edge signs, kink numbers and trapped areas 
are indeed topological invariants (ie, they are invariant under continuous de- 
formations of n that preserve the tangent boundary conditions) and that two 
tangent unit-vector fields on P are homotopic if and only if their invariants 
are the same. For convenience we have slightly adapted the notation of |1] 
to suit the case of prisms (the conventions, however, are the same). 

In what follows we restrict our attention to a subset of the allowed prism 
topologies which we call reflection symmetric. Let v and w denote a pair of 
vertices related by a refiection though a midplane of the prism (and therefore 
joined by an edge). For refiection-symmetric topologies, the edge signs at v 
and w are the same while the kink numbers and trapped areas differ by a 
sign. That is, 

e] = ej, kJ = -kY, n^ = -n'^. (8) 

It follows that at vertices related by two reflections (ie, at diagonally opposite 
corners of a face), the invariants are the same, while the invariants of vertices 
related by three reflections (at diagonally opposite corners of the prism) are 
related as in (jHl). 

For reflection-symmetric topologies the invariants are determined by their 
values at a single vertex. For deflniteness we take this vertex to be the 
origin, and henceforth denote the invariants simply as {e,k,Q). The surface 
separating the origin from the other vertices will be denoted by C. It is 
straightforward to check that (jH)) implies that the sum rules Q and ((Zj) are 
automatically satisfled. 
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The terminology stems from the fact that every reflection-symmetric 
topology has a reflection-symmetric representative, ie a configuration n which 
is symmetric under reflections through the mid planes, 

n(x, z) = n(L^ - x, y, z) = n(x, Ly -y,z) = n(x, y, L^- z). (9) 

Let 

i? = {r|0<r, <iL,} (10) 

denote the octant of the prism with the origin as vertex. Then a reflection- 
symmetric configuration is determined by its values in R. It is straightfor- 
ward to verify that (jH)) implies the relations (jH)). Conversely, given a configu- 
ration n' with reflection-symmetric topology (e, /c, f2) but which is not itself 
reflection symmetric, we can construct a reflection-symmetric configuration 
n satisfying (jH)) with invariants (e, fc, Vt) (just take n = n' in the prism octant 
R and define n elsewhere via ©)• 

In we introduced certain additional integer-valued topological invari- 
ants, called wrapping numbers. As the preceding discussion implies, the 
wrapping numbers are not independent of the edge signs, kink numbers and 
trapped areas, but rather can be expressed in terms of them. We briefly 
recall the definition and properties of the wrapping numbers, as they are 
central to the discussion to follow. 

Let s denote a regular value of n restricted to C. That is, on C, there is 
a finite number of points where n takes the value s, and, at any such point, 
the Jacobian of the map nc : C — >^ S*^ is nonsingular, so that, from (jH)), 
(D ■ C)(s) 7^ 0. The wrapping number at s, denoted w(s), is a signed count 
of the pre-images of s on C, denoted r^, the sign determined by whether ric 
is orientation-preserving (+) or reversing (— ) at r^. Thus 



5^sgn[(D.C)(r,)]. (11) 



To express the wrapping number in terms of the other invariants jH E] , let 
Upe denote the disk of radius e about on C, and let C — Up^ denote 
C with these disks excised. Let dC denote the boundary of C and dUp^ the 
boundary of Up^. Choosing s as the south pole of the polar angles (6*, 0) in 
(jHI) and using Stokes' theorem, we get that 

Q = \im sin 6 d6 A d(f) = i - lim / )(l-cos 



'^{C-ZpUp,) \Jn{dC) '^^Jn^ZpdUp 



(12) 

T^i^pdUpe) consists of p small circuits about s, and the integral over these 
circuits in p2|l gives, in the limit e — > 0, 47r times w{s). n{dC) consists 
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of the spherical right triangle with vertices e^j, j = x, y, z, along with kj 
circuits of the great circle normal to j. Each great circle contributes ±2tt to 
the integral in ()12|1 according to its orientation with respect to the polar axis 
through s, while the spherical triangle contributes ±7r/2 (its signed area) or 
±77r/2 according to whether or not it encloses s. Keeping track of signs one 
gets 

1 ^ 1 , I —I, if o'j = sgne,- for all ?, , , 



otherwise. 



where aj = sgnsj. From ()13|1 it is clear that w{s) is a topological invariant. 
In fact, w{s) can be defined so long as s does not lie in a coordinate plane 
(ie, even if s is not a regular value of n) as the degree of a certain continuous 
^2 _^ g2 j-Qg^p constructed by gluing n : C — S"^ to a reference map which 
coincides with n on the boundary dC 

(fT^ also implies that w{s) depends only on the signs of the components 
of s, ie on the octant of 5*^ to which s belongs. In what follows, we label 
octants by a triple of signs a = {crx,cry,az), so that denotes the octant 
{s|sgnsj = CTj}. For convenience we let w„ denote the value of w{s) for 
s G Oa- Summing over octants in (fT^ . we get that 

^=I^5Z^- (14) 

(T 

(the terms in (jlHj) involving ej and kj cancel in the sum). 



3 Lower bound for the elastic energy 

In the continuum theory of nematic liquid crystals [12], the elastic, or Frank- 
Oseen, energy of a configuration n is given by 

E{n) = j [i^'i(divn)^ + ii'2(n ■ curl n)^ + 7^3 (n X curl n)^ 

+ /s:4div((n- V)n- (divn)n)] (15) 

Tangent boundary conditions imply that the contribution from the i^4-term, 
which is a pure divergence, vanishes. In the so-called one-constant approxi- 
mation, the remaining elastic constants Ki, K2 and are taken to be the 
same. In this case, (fT3j) simplifies to 
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We shall use the one-constant approximation in what follows. 

Let E^i^{e, k, fl) denote the minimum (infimum) energy for configurations 
with reflection-symmetric topology {e,k,Q). In jZj we obtained the lower 
bound 

^min(e,A;,^]) >8L,|^]|. (17) 
In view of p4p . this may be written as 



£'min(e, k, Q) > AtiL^ 



Here we derive a new lower bound which, in general, is an improvement on 

m- 

Theorem 3.1. 

{e,k,VL)>AnL^^\w^\. (19) 

Proof. Let n be a configuration with reflection-symmetric topology (e, fc, VL) 
for which the energy p6|) is finite. As shown in FB|, we can, without loss 
of generality, take n to be smooth (smooth configurations are dense in the 
space of finite-energy configurations with respect to the energy norm). 

We can assume that the energy of n in i? is not more than its energy in 
any other octant of the prism (we can replace n(r) by n(7^ ■ r) for a product 
TZ of reflections through mid planes; the reflected configurations have the 
same topology and energy as n). Then 

E(n) > 8 /" iVnprfV > 8 /" |VnprfV, (20) 

jR Jr<L^/2 

where the last integral is taken over the positive octant of the ball of radius 
Lz/2 about the origin. Using the local inequality for the energy density 

IHIEIIZI, 

(Vn)^ > 2|D| > 2|D • f|, (21) 

we get that 

r rLz/2 r 

E(n) > 16 / |D(r) ■ f|rfy = 16 / dr |D(r) • f | (22) 

Jr<L^/2 Jo JreCr 

Here, Cr is the positive octant of the sphere of radius r about the origin, 
with area element dSr- 
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We partition Cr into preimages of the octants O^- of 5^, writing 

E(n) > 16 f ' dr f ( V f ds (552(n(r), s) | |D(r) • f | rfS,. (23) 

Here ^52 is the normahsed Dirac deha-function on S*^, so that ds 632 (n(r), s) 
equals one if n(r) G 0^ and is zero otherwise. We interchange the integrals 
over r G C,. and s and take the absolute value outside these integrals to 
obtain 



E(n) > 16 



L,/2 



ds 



552(n(r),s)D(r) • f 



reCr 



(24) 



For s a regular value of n (by Sard's theorem, regular values are of full 
measure), we get that 

/ <552(n(r),s) dSr = J] \det dncA^pT' = J] |D(rp) ■ ff^ , (25) 



where the sum is taken over the preimages G C,. of s, and we have used 
(0)- Substituting into (j^ . we get that 



(26) 



E{n)>16 dry2 ds V sgn (D(rp) ■ f ) 

-'0 „ Jo„ 



From (fTT|) . the sum over p is just the wrapping number w{s) = w„, so that 
the integral over s trivially gives a factor of 71/2 (the area of Oo-). Then 
the integral over r trivially gives a factor of Lz/2. The required result (fTIHl 
follows. 

□ 



4 Conformal and anticonformal topologies 

The new bound (fT^ agrees with the previous bound (fTHj) for topologies 
where Yla ~ Eo-'^^crl' where the nonzero wrapping numbers all have 
the same sign. We will say that a reflection-symmetric topology is conformal 
if Wa ^ for all a, anticonformal if w„ > for all a, and nonconformal 
if neither of these conditions holds. Thus, the new bound constitutes an 
improvement for nonconformal topologies. 

It is useful to characterise the conformal and anticonformal topologies 
directly in terms of the invariants (e, k, Q). 
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Proposition 4.1. Define functions fi;^(e, k), where x = ±; follows: 
For x^xeyC^ = 1, 



7 



^ / , N , I -t-7: if K^ikj < for all ?, , 

fi^(e,/c) = 27rj]|/c,| + 27r<^ f (27) 
■'^ — ' ' — otherwise. 



For xe^e-ytz 



4' 
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^ / , N , I -i-7, if ye-ikj < for all ?', 

Q^{e,k) = 27rJ2\k,\-2n\ f ^ ^ ^, (28) 
— ' ' — otherwise. 



4' 



r/ien i/ie reflection- symmetric topology (e, fc, fi) zs conformal if and only if 
Q < — r2_(e, fc) and anticonformal if and only if ^ > Vt^{e.,k). If equality 
obtains, ie Q = ~Q^{e,k) or Q = Q^{e,k), then at least one wrapping 
number must vanish. 

Proof. The condition xu!a > for all a is equivalent to (e, k, Q) being con- 
formal (x = — ) or anticonformal (x = +)• From (jlHj) . > for all a if 
and only if x^ > ^x(^' where 



^^^(e, k) = 27rmax | 



-|, if cTj = Cj for all j 
-i, otherwise 



(29) 

with x^ = ^x(^' ^) ^'^'^ only if = for some a. In fl29|l . to realise 
the maximum we may take, for all j such that kj ^ 0, aj = — xsgn/cj, and 
thereby replace —x<^jkj by \kj\ for all j. Thus, 

^ . , . . , . I +l7r, if o",- = Cj for all j, . . 

\L^[e,k) = 2ti y \kj\+ max xe^ej^e^ x < ^ j ■> _ (^39 j 

"~ '^j|fcj=o [^l^' otherwise 

It remains to maximise the second term in (j30|) with respect to the (Xj's 
for which kj = 0. Suppose that x^x^yCz = 1. Then, provided = cj for 
all kj 7^ 0, ie provided — xsgnfcj = cj for all kj 7^ 0, the maximum value 
attained by the second term in (jHUj) is 7tt/2. Otherwise, the maximum is 
— 7r/2. This is in accord with (j^Tj). Next, suppose that x^x^y^z = —1- The 
maximum value attained by the second term in (j3fl|) is 7r/2 unless all the fcj's 
are nonzero and — xsgnfcj = cj, in which case the maximum is —7tt/2. This 
is in accord with ()28j) . □ 
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In [7j we introduced certain reflection-symmetric configurations in P 
which we called conformal and anticonformal. We show next that the confor- 
mal and anticonformal topologies are precisely those which have conformal 
and anticonformal representatives. To proceed, we briefly recall the proper- 
ties of conformal configurations (fS^) (as discussed below, the treatment of 
anticonformal configurations is analogous). A reflection-symmetric configu- 
ration n is said to be conformal if, in the prism octant R, i) n is radially 
constant, ie n(Ar) = n(r), and ii) n is conformal, ie the map t Vjn(r) 
from vectors t orthogonal to f to vectors Vtn(r) orthogonal to n(r) preserves 
orientation, angles and ratios of lengths (or else vanishes). 

Conformal configurations are conveniently represented via stereographic 
projection as analytic functions f{w), 

The domain of fiw) is the quarter-unit-disk Q given by |w| < 1, < Kew < 
1 and < Imw < 1. The boundary of Q consists of the real interval 
< w < 1 (which corresponds to the xz-i&ce of i?), the imaginary interval 
< —iw < 1 (which corresponds the yz-i&ce), and the circular arc |w| = 
1, where < axgw < 7r/2 (which corresponds to the xy-face). Tangent 
boundary conditions imply that i) f{w) is real for w real, ii) f{w) is imaginary 
for w imaginary, and iii) |/(w)| = 1 if \w\ = 1. Assuming that f{w) has 
a meromorphic extension to the extended complex plane, these conditions 
imply that if is a zero of /, then —w^, and w^, are zeros, while is a 
pole. The meromorphic functions which satisfy these conditions are rational 
functions of the following form: 

j=l ^ J ' k=\ ^ ' 

Here, e = ±1 and n, an odd integer, gives the order of the zero or pole 
of / at the origin, a is the number of zeros of / (p-,- = 1) and poles of 
/ (pj = —1) on the real interval (0,1), with positions ordered so that 
< Ti < ■ ■ ■ < Ta < 1. Similarly, h is the number of zeros of / (ufc = 1) 
and poles of / (ajt = —1) on the imaginary interval (0,i), with positions is^ 
ordered so that < Si < ■ ■ ■ < < 1. Finally, c is the number of zeros of / 
(t; = 1) and poles of / [ji = —1) in the interior of Q, with positions ti. 
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The edge signs, kink numbers and trapped area of conformal configura- 
tions are given by 



= e(-l)^ Cy = e(-l)^(-l)("-i)/2^ = sgn n, (33) 
ky = -|(-l)°e. ^ " (-1)'')^^) ' (34) 

a be 

j=\ k=i 1=1 



and 



= +2(a + 6) +4c)7r. (35) 

As explained in the expressions for the edge signs and the trapped area 
are easily derived (in particular, ()35|) follows from consideration of the degree 
of / on the extended complex plane). A derivation of the expression for the 
kink numbers, which was deferred in [7|, is given here in Appendix 1X1 

Clearly, a conformal configuration has a conformal topology; the orientation- 
reversing property ensures that the wrapping numbers cannot be positive. 
Below we establish the converse fact; every conformal topology (e, k, Q) has 
a conformal representative. The demonstration splits into four cases accord- 
ing to the sign of e^eyCz and of ejkj. In each case we exhibit the parameter 
values, expressed in terms of e, k and Q, for a particular conformal configura- 
tion. It is then straightforward to verify - we omit the explicit demonstration 
- that the specified parameters are admissable (ie, that n is an odd integer; 
a, b, c are nonnegative integers; e, pj, ak, ti are signs), and that, with these 
parameters, the values of the invariants given by (jHHj) - (jH^j) are just (e, k, Q). 
Deriving the exhibited values involves a systematic and slightly tedious inves- 
tigation of (jSni)^ For the sake of brevity, these details are also omitted. 

Case la. e^eyCz = 1 and ejkj > for all j . Let 

a = 2\ky\ - 1, pj = {-lycz, 
b = 2\kJ-l, ak = {-\fez 




-^fi- 1^.1-1^,1 + 1, n = {^-- l^'^^^l' (36) 



/ > \kzy 
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Note that (fT^ implies that c — (|/c^| — 1) is nonnegative and even. Here and 
in the cases to follow, we do not specify the positions of the zeros and poles 
explicitly. 

Case lb. exCyCz = 1 and Cjkj < for some j . Without loss of generality, we 
may assume that j = z. This follows from considering the fractional linear 
transformation 

r(w) = "—^ (37) 

which maps Q onto itself while cyclically permuting its vertices (r corre- 
sponds to the 27r/3- rotation on S'^ about the axis (1,1,1)). Therefore, if / is 
a conformal configuration, so is / given by 

~f = rofor-\ (38) 

It is easily verified that / and / have the same trapped areas while their edge 
signs and kink numbers are related by cyclic permutation. 



), h — (/c^, ky). (39) 



6 \(^zi ^xi (^y 

Letting e^fc^ < 0, we take 

e = Cx, n = (4|A;^| + 1)6^, 
a = 2\ky\, pj = -{-lye^sgnky, 
b = 2\kx\, o-fc = -(-l)^ej^sgn /c^, 

c=-^n-\k.,\-\ky\-\k,\-l, n = {-iy. (40) 

Note that ()13j) implies that c is nonnegative and even. 

Case 2a. e^eyCz = —1 and Cjkj < for some j . As in Case lb, without loss 
of generality, we may take e^k^ < 0. Let 

€ = 6^, n= -{Akz + Cz), 
a = 2\ky\, Pj = -{-lye^sgnky, 
b = 2\ka:\, ak = -{-!)'' CySgnkx, 
c = - \k.\ - \ky\ - \kz\ + i, n = i-iy. (41) 

Note that (fT^ implies that c is nonnegative and even. 
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Case 2b. e^eyCz = —1 and ejkj > for all j. Let 

a = 2\kyl pj = e^i-iy 
b = 2\k^\, ak = e^{-l 



Note that (jlHj) implies that c — (I/e^I + 1) is nonnegative and even. 

Anticonformal configurations are given by antianalytic functions, in anal- 
ogy with the conformal case. Specifically, if f{w) is a conformal configuration, 
then f{w) is anticonformal with invariants (e, k, Q) given by 

6 = (cj;, Cy, e^), k = ( /cj:, /Cy, kz), ^ = fi. (43) 



Thus, conformal and anticonformal configurations are in one-to-one corre- 
spondence. Also, given any {e,k,fl) and {e,k,fl) related as in PHj). one can 
verify from (j2ZI) and ^ that Q+{e,k) = Q-(e,k), so that Q > Q+{e,k) if 
and only if < — fi_(e, fc). Thus, conformal and anticonformal topologies 
are in one-to-one correspondence, and representatives of every anticonfor- 
mal topology may be obtained from complex conjugation of the associated 
conformal representative. 

The preceding discussion may be summarised as follows: 

Proposition 4.2. A reflection- symmetric topology is conformal if and only 
if it contains a conformal configuration, and anticonformal if and only if it 
contains a anticonformal configuration. 



5 Upper bound for elastic energy 

In j7] we showed that for a conformal or anticonformal topology (e, k, Q), 

E^i^{e,k,Q) <8L\Q\. (44) 

where 

L = {Ll + Ll + Ll)y^. (45) 

We note that the upper bound differs from the lower bound ()17p by a factor, 
L/Lz, which depends only on the aspect ratios of the prism, and not on 
(e, k, fl). 
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Here we derive an analogous upper bound for nonconformal topologies. 
To this end, we construct representatives n. As in the conformal and anticon- 
formal cases, we take these to be reflection-symmetric and radially constant 
in R. In i?, n is taken to be of the form 



^ ~'~ ^ ' ''r) = F{w,w), where w = — (46) 



1 + riz J r + z 

We take F to be a juxtaposition of analytic and antianalytic domains (in 
which the local estimate ()21|) for the energy density becomes an equality), 
separated by an interpolating domain of small energy. For definiteness, we 
take n < (the case f2 > is treated analogously). Let / denote the confor- 
mal configuration with topology (e, k, — f2_(e, k)), so that / is the conformal 
configuration with the largest trapped area compatible with e and k. Let wq 
denote a point in the interior of Q, and let D^{wo) = {w\\w — Wo\ < e} de- 
note the open e-disk about Wq. Choose Wq and e so that D2e{wo) is contained 
in Q and contains no poles of /. Let 

W=l{n + n.ie,k)). (47) 

Since (e, k, Q) is nonconformal, is a positive integer. We let 

{/(w), weQ-D2eiwo), 
sfiw) + (1 - s)if{wo) + {w- wo)"^), w E D2e{wo) - D^iwo), 
f{wo) + e'^^iw - wo)-"^, w e D,{wo), 

(48) 

where 

s[w,w) = (49) 

e 

(so that s varies between and 1 as |w — Wo| varies between e and 2e). 
Thus, in Q — D2e{wQ), F coincides with / and therefore is conformal, while 
in Df{wo) it is anticonformal. 

Let us verify that F has the required topology. Since it coincides with / on 
the boundary of Q, F has the same edge signs and kink numbers as /, namely 
e and k. As for the trapped area, from ^ and (jlUI) it is straightforward to 
derive the general expression 

^(^)=/^4 '^"^^'^yj^-f'' ^V (50) 

()50|1 can be evaluated by dividing the domain of integration as in fl48j) . The 
contribution from Q — D2e{wQ) is, to O(e^), just the trapped area of /, namely 
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—Q-{e,k) (substituting / for F in (j^Uj) . the contribution from Z}2e(^o) is 
0(e^)). Consider next the contribution from the disk D^{wo). Here F = 
/(wo) + e"^^ {w — wo)~^ , so that F covers the extended complex plane, apart 
from an e'^-disk about f{wo), W times with positive orientation. It follows 
that the contribution to (jSOI) is, to within 0{e^) corrections, inW. The 
remaining contribution, from the annulus D2e{wo) — D^{wo), is O(e^). This 
is because the area of the annulus is O(e^), while the integrand in (j3n|) may 
be bounded independently of e (by assumption, / has no poles in D2e{wo)). 
Since the trapped area is an odd multiple of 7r/2, it follows that, for small 
enough e, 

n{F) = -n^{e,k) + 4:7iw = n. (51) 

By estimating the energy of the nonconformal representatives we can 
obtain the following upper bound for £^min(e, k, Q): 

Theorem 5.1. Let {e,k,Q) denote a nonconformal topology. Then 

{e,k,VL) <?,QtiL^\w^\. (52) 

Proof. From (fT^ and ()46|) one can derive the general expression for the 
energy, 

E{F) = 16 [ 4|r(^)| '^"f''|j^:f'' dV (53) 



where r{w) is the point on the boundary of R which has w as its stereographic 
projection. Since |r(t(;)| < L/2, it follows that 

E{F) < 8LA{F), (54) 

where 

■4(F)^/4 l"fi:+lf:f^ V (55) 



A{F) represents the unoriented area of n{Cr). The expression for A{F) 
differs from the expression (jHn|) for Q{F) in the relative sign of the w- and 
w-derivative terms. (Thus, for conformal and anticonformal configurations, 
one obtains the estimate ^^■) Arguing as for (jSTj), we have that 

A{F) <\n4e,k)\+A7rW. (56) 

From ()14|1 . this may be written as 

A{F) <AttW -\7iJ2^<r-, (57) 
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where i^o-- are the (nonpositive) wrapping numbers of /. From (fT^ and (^Tj), 
the W(j-'s are related to the wrapping numbers of F, denoted according 
to 



Substituting into (jHTjl . we get that 

cr a 

One easily establishes the estimate 

Indeed, since f2 < by assumption, it follows from (fT^ that 

< 0. 



(5^ 



(59) 



(60) 



(61) 



Since / has trapped area —Q^{e,k), it follows from Proposition 14.11 that 
there is at least one octant, say cxo, in which / has zero wrapping number. 
From (jKH|) . Wa-o = W. Let Yl'a denote the sum over octants with ctq omitted. 
Then 

Y^w^ < -W. (62) 

cr 

It follows that 

= Y^M + W > + W>2W. (63) 

(T cr cr 

Substituting into (fKTjl . we get that 

A{F) <9x ^ttJ2M- (64) 

cr 

The required result, follows from substitution into (jSH)- □ 

6 Discussion 

For nonconformal topologies, the ratio of the upper and lower bounds for 
Erain, as given by and (fT^. is dL/L^. By finding representatives of 
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lower energy, it might be possible to obtain a ratio closer to the confor- 
mal/anticonformal result, L/L^ (which can be further improved by more 
accurate estimates of the energy of the representatives [Zj). 

In IHj, we described, for conformal and anticonformal topologies, a transi- 
tion in topologically nontrivial equilibrium (infimum-energy) configurations, 
from singular, in the case of a cubic domain, to smooth, as the prism as- 
pect ratios are varied. Singular configurations, when they appear, are limits 
of configurations which differ from the topologically simplest "unwrapped" 
configurations in thin tubes along the prism edges. It would be interesting 
to investigate whether similar transitions occur for nonconformal topologies. 
The nonconformal representatives differ from conformal/ anticonformal con- 
figurations only in a tube (a disk in the two-dimensional stereographic de- 
scription (HH))), and depending on the aspect ratios, it may be energetically 
advantageous for these tubes to collapse to edge singularities, or not. 

AM was supported by an EPSRC/Hewlett-Packard Industrial CASE Stu- 
dentship. MZ was partially supported by a grant from the Nuffield Founda- 
tion. We thank CJ Newton and A Geisow for stimulating our interest in this 
area. 

A Kink numbers of conformal configurations 

Taking n to be a conformal configuration with stereographic projection / 
given by ()32j) . we derive formulas for the kink numbers k = {kx,ky,kz) in 
terms of the parameters of /. 

Formulas for k^ and ky. For definiteness, consider the calculation of ky. Let 
rl^), < r < 1, denote a small quarter-circular arc on the xz-face of R 
starting on the z-axis and ending on the x-axis (so that r(r) is positively 
oriented with respect to the outward normal — y through the centre of the 
face). Let ^'(r) = n(r(r)). Then ^'(t) describes a curve on S"^ along the 
great circle in the a;2;-plane, starting from and ending at e^^x. ky is the 
winding number of ^'(t) relative to the shortest arc joining e^z to e^jX, with 
anticlockwise windings about — y taken as positive, ky is given by the number 
of times ^'(t) crosses a given point, say z, counted with a sign according to 
orientation. Let Tp denote the parameter values at these crossings. Assuming 
that v'ijp) 7^ 0, we get that 

K = -Y. (^'(^p) • + \{l + e,) ■ i(e, - sgn {u'{Q) ■ x)). (65) 

Tp>0 

If = 1 then i/(0) = z; the last term in accounts for a possible contri- 
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bution in this case. There is no contribution if sgn (^''(0) ■ x) = Cx-, as for the 
shortest arc joining z to e^x (for which fc^ = 0). Otherwise, the initial point 
constitutes a crossing with sign e^. 

Under stereographic projection, r(r) corresponds to w{t) = r, the cross- 
ing vi^Tp) = z corresponds to /(rp) = 0, and sgn {i^'{Tp) ■ x) corresponds to 
sgn f'{Tp). The zeros of / on the real interval (0, 1) are given by the r^-'s with 
Pj = 1. Thus, becomes 



J2 sgn/'(r,) + 1(1 + e,) ■ i(e. - sgn/'(0)). 



(66) 



i|pj=i 



For simplicity, let us assume that the r^-'s are all distinct and ordered so that 
< ri < • ■ ■ < < 1. In this case, for pj = 1, f'{rj) ^ 0, and from (jS21), we 
have that 



sgn/'(rj) = sgn 



2r. 



in 



2 _ 2 \ P" 



1 \'m'j 



ef-lV. 



(67) 



If = 1, n is positive, and ti? = is also a zero of /. As /'(O) vanishes if 
n > 1, we replace sgn/'(0) by lim^^o sgn /(w) = e. (pT?)|) becomes 



(6^ 



Recalling that = e(— 1)", with some further straightforward manipulation 
we obtain 



(69) 



which is just the expression given in ()H4j) . In fact, ()69j) holds even if some of 
the r/s coincide. 

The expression for is similarly derived, with r(r) taken to be a quarter- 
circular arc on the zy-ia.ce of R with projection w = i{l — t). Details are 
omitted. 



Formula for k.. Let r(r), < r < 1, denote a small quarter-circular arc 
on the xy-face of R starting on the x-axis and ending on the y-axis (so that 
r(r) is positively oriented with respect to the outward normal — z through 
the centre of the face). Under stereographic projection, r(r) corresponds 
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to w{t) = exp(i7rr/2), and n(r(r)) to /(exp(z7rr/2)). is the winding 
number of /(exp(z7rr/2)) on the unit circle in the complex plane relative to 
the shortest arc joining to Cyy. Clockwise windings are taken as positive. 
It follows that 

kz = (arg/(exp(m/2)) - arg/(0)) + ^ sgn (e^.ey), (70) 

where arg /(exp(i7rr/2)) is taken to be continuous in r and the last term 
ensures the winding number is zero for the shortest arc joining to Cyy. 
Referring to (jS2I), each factor of the form [{w"^ / {p^w"^ ± 1)]^ contributes 
,^7r to the change in arg/ in (f7n|) . while contributes mT/2. Thus we get 

a be 

kz = -\^p3-\^(rk-^ri- \{n-exey), (71) 
j=i k=i 1=1 

as in dSl). 
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